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(Dated: February 3, 2008) 

We present the results of a search for the decay 5° —* t + t~ in a data sample of (232 ± 3) x 10 6 
T(4S) — * BB decays using the BABAR detector. Certain extensions of the Standard Model predict 
measurable levels of this otherwise rare decay. We reconstruct fully one neutral B meson and seek 
evidence for the signal decay in the rest of the event. We find no evidence for signal events and 
obtain B(B° -> t+t~) < 4.1 x 10" 3 at the 90% confidence level. 

PACS numbers: 13.20.He, 14.40.Nd, 14.60.Fg 



None of the leptonic decays B° — > (£ = e, /x, t) 

have been observed. In the standard model of particle 
physics, the decays can be mediated by box and penguin 
diagrams (Fig.|IJ. The standard model produces only the 
combinations £^i£ and t\£~fi- The amplitudes for the 
decay of a spin-zero particle to these states are propor- 
tional to mi and thus the decay rates are suppressed by 
(mi/ms) 2 . The suppression is smallest for B° — > t + t - 
due to the large r mass. The standard model prediction 
for the B° — > t + t~ branching fraction is 1] 



W 



W+ 



1't 




FIG. 1: Standard model box and penguin processes that can 
mediate B° -> £ + T (q = t, c, u). 
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where fs is the B decay constant and Vtd is the Cabibbo- 
Kobayashi-Maskawa matrix element. The theoretical un- 
certainty on fg and the experimental error on Vtd domi- 
nate the uncertainty on the predicted branching fraction. 

Extensions of the standard model containing lepto- 
quarks, which couple leptons to quarks, predict enhance- 
ments for B(B° -> T+T-) that are proportional to the 
square of the leptoquark coupling. In theories that con- 
tain two Higgs doublet fields, the rate can be enhanced by 
powers of tan /3, the ratio of vacuum expectation values of 
the two Higgs doublet fields [10]. Since B° -> £+£~ has 
not been observed, one can only constrain model parame- 
ters using the measured branching fraction limits. While 
tan/3 is constrained by all three modes (£ — e, fj,, t), only 
B° — > t + t~ can constrain the coupling of a leptoquark 
to the third lepton generation or other new physics in- 
volving only the third generation. 

The analysis described here provides the first upper 
limit on £>(B° — > t + t~). The data were collected with 
the BABAR detector at the asymmetric PEP-II e + e~ stor- 
age ring. A full description of the BABAR detector is given 
in Ref. 5]. In brief, charged-particle momenta are mea- 
sured with a tracking system comprising a silicon ver- 
tex detector (SVT) and a drift chamber (DCH) placed 
within a highly uniform 1.5-T magnetic field generated 
by a superconducting solenoid. Electron and photon en- 
ergies are measured with an electromagnetic calorimeter 
(EMC) constructed with Thalium-doped Csl scintillating 



crystals. Muons are distinguished from hadrons in a steel 
magnetic-flux return instrumented with resistive plate 
chambers (IFR). Charged particle identification is pro- 
vided by a Cherenkov detector (DIRC) and the tracking 
system. The data sample consists of 210 fb _1 collected at 
the peak of the T(45) resonance, which corresponds to 
232 ± 3 million BB pairs. The expected background and 
the expected signal efficiency are obtained from Monte 
Carlo simulation samples. The sample events were gen- 
erated with the EvtGen event simulator 6} and propa- 
gated through a detailed model of the BABAR detector 
using the GEANT4 detector simulator 0. 

Isolating B° — > t + t~ poses a unique challenge. This 
decay contains at least two and as many as four neutrinos, 
so there is no kinematic discriminant that separates sig- 
nal from background due to undetected particles. Since 
two B mesons are produced in an T(AS) decay, the mis- 
assignment of decay products to the parent B must be 
avoided. We completely reconstruct one B candidate in 
each event (hereafter referred to as the companion B) 
and search for the signal decay among the remaining de- 
tected particles. The combinatorial background in the 
companion-B reconstruction is determined by a fit to 
the companion-B invariant mass distribution. We em- 
ploy the parameters 

™ ES = ^E* b l im -p*£ (2) 
AE = E% - El eam , (3) 

where p* B and E B are the reconstructed companion-B 
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momentum and energy in the center-of-mass (CM) frame. 
E£ eam is the beam energy in the CM frame. The m ES 
distributions are fit with a probability density function 
composed of a Crystal Ball function || to model the peak 
at the B mass and an ARGUS function to model the 
nonpeaking combinatorial background. 

The companion B is fully reconstructed in a hadronic 
mode B° -> D&X , where flW is either a D+ 0, -D , 

ting < 

1 1 1| . D* + mesons are re- 
constructed in the channel D°ir + . D° mesons are recon- 
structed in the channels K~tt + , K~tt + it , K~ir + ir~ir + , 



or D* + and X is a system consis ting of up to five particles 
of the type ir ± , tt°, K ± , or K%' 



D + mesons are reconstructed in the S 



and K s ir + ir~ 

channels K s tt + , K~tt + tt + , K s tt + tt°, K s w + ir + ir~ , and 
K + K~ir + . The AE of the companion B is required to 
be within two mode-dependent standard deviations of 
the mean when no ir° is present, or to satisfy —0.09 < 
AE < 0.06 GeV for reconstructions with one or more 
7T°. If more than one B candidate is reconstructed in the 
same mode, the reconstructed B with the smallest \AE\ 
is selected. For each mode, the purity B pur is the ra- 
tio of the number of events before signal selection in the 
fitted peak to the total number of events in the region 
5.27 < m ES < 5.29 GeV. Only events reconstructed in 
a mode with B pur > 0.12 are selected, which results in 
the reconstruction of 147 distinct modes in the data sam- 
ple. If B candidates are reconstructed in more than one 
mode, the B reconstructed in the mode with the highest 
B pur is selected as the companion B. 

We estimate the total companion- B yield from all re- 
constructed modes using the BB and qq (q — u, d, s, c) 
simulated samples before applying the signal B° — ► t + t~ 
selection. We first remove the peak from the B°B° 
simulated sample using the fitted Crystal Ball probabil- 
ity density function. Subtracting the simulated combi- 
natorial background m ES shape, fitted to the data be- 
low 5.26 GeV, from the data distribution yields a nomi- 
nal companion-B yield of N B a B o = (2.80 ± 0.27) x 10 5 
(Fig. |2| . The systematic error on N B o B o is estimated 
to be 10% by varying the fit region and by varying the 
combinatorial background composition with event-shape- 
variable cuts. 

The companion-B decay products are removed from 
the event and the signal-i? characteristics are sought 
among the remaining particles. The dominant back- 
ground to B° — > t + t~ arises from decays b — > W~c(—* 
W + s), in which the s quark hadronizes into a K L that 
escapes detection and the virtual W + and W~ mimic 
the virtual W + and W~ emitted by the signal r. A sec- 
ondary background originates in events in which two op- 
positely charged particles are lost outside the detector 
fiducial region. We select signal events that are con- 
sistent with each r decaying to a single charged par- 
ticle (and one or two v) by selecting events with zero 
net charge and two tracks in the recoil system. Each 
track must leave at least twelve hits in the DCH, origi- 
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FIG. 2: Above, the m B s distribution for the hadronic com- 
panion B in data (dots) and scaled simulated background 
(upper histogram) before the signal B° — > r + r~ selection 
is applied; the lower histogram is obtained by subtracting 
the background from the data. The companion- B yield is 
N b obo = (2.80 ± 0.27) x 10 5 . Below, the m ES distribution 
after the signal B° — » t + t~ selection. The fitted proba- 
bility density function (short-dash) and its ARGUS compo- 
nent (dash) are superimposed on the data (dots). We obtain 
N b s = 263 ± 19 events in the peak. 



nate within 10 cm of the beamspot in the beam direction 
and within 1.5 cm in the transverse direction, and have 
a transverse momentum of at least 0.1 GeV. To elimi- 
nate background originating from b — ► W~c(-^ W + s) 
events, the selection rejects events with identified K + , 
K° s , or K° L . The K± candidates are identified by a neu- 
ral network with inputs taken from the SVT, the DCH, 
and the DIRC. The Kg candidates are identified as a 
7r + 7r~ pair with invariant mass consistent with the Kg 
mass (0.473 < m v+7T - < 0.523 GeV). The K° L candidates 
are identified from clusters in the EMC that have not 
been associated with a charged track or included in a 
candidate 7r°. A neural network is employed to identify 
K\ candidates using the cluster energy and shower-shape 
variables, which discriminate hadronic from electromag- 
netic showers. 

The multiplicities of e, [i, and ir° in the recoil sys- 
tem must be consistent with each r decaying in one of 
the channels r — > irv,pv, evv or pvv (Table P). The e 
candidates are identified with dE/dx measurements from 
the DCH and shower shape variables from the EMC. 
The y, candidates are identified with variables from the 
IFR (to reject the ir hypothesis) and EMC (to reject 
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TABLE I: Signal B 1 — » r + r branching fraction and require- 
ments by mode [l — e,p). 



Selection Mode 


B(%) [12] JVe + iV M 


2V 


m n7T o 


T + T~ — + lvv/l'vv 


12.4 


2 





NA 


T + T~ — > ivv j-KV 


7.8 


1 





NA 


T + T~ — > IVV / pV 


17.7 


1 


1 


[0.6, 1.0] GeV 


t + t~ — > nv/nv 


1.2 








NA 


t + t~ — » nv/ pv 


5.6 





1 


[0.6, 1.0] GeV 


T + T~ — » pV 1 pV 


6.3 





2 


[0.6, 1.0] GeV 



TABLE II: e S i 3 , N expecte d and N b s obtained from individ- 
ual fits by signal mode. The errors are statistical and fit 
error added in quadrature. Branching fractions are included 
in the efficiency estimates. The irv/nv channel is dominated 
by crossfeed from other signal channels. 



Selection Mode 


e S i S (%) 


N expected 


N obs 


t + t~ — > lvv/l'vv 


0.9 ±0.2 


46 ±4 


54 ± 7 


T + T~ — ► Ivv/lXV 


1.5 ±0.3 


122 ± 6 


105 ± 11 


T + T~ — > TTv/lVV 


1.5 ±0.3 


89 ±6 


80 ±11 


T + T~ — > pV 1 pV 


0.3 ±0.1 


21 ±3 


15 ±6 



the e hypothesis). Track candidates that are not iden- 
tified as e, /j or if are assumed to be tt. Events with 
7T° are vetoed unless the n° can be associated to a n + 
such that the invariant mass is consistent with the p 
mass (0.6 < m v + v o < 1.0 GeV). The 7r° candidates are 
formed from pairs of 7 candidates with invariant mass 
0.090 < to 77 < 0.170 GeV, with each 7 having an energy 
greater than 0.030 GeV. Since the presence of residual 
unassociated energy in the EMC (E res ) is a strong indi- 
cation that an unreconstructed ir° or K° is present, we 
require E res < 0.11 GeV. 

The r-daughter candidates are Lorentz-boosted with 
the companion-_B momentum. While distributions of 
the momenta p + and p_ of the charged daughters ex- 
hibit no discrimination from the background momen- 
tum distributions, correlations among |p+|, |p-|, and 
cos# = p+ ■ p_/|p+||p_| afford some discrimination, es- 
pecially when categorized by signal B° — > t + t~ selec- 
tion mode. Cascade decay background events manifest 
an asymmetry in |p+| and |p_| that is not present in sig- 
nal events. The parameters |p+|, |p-|, cos#, E res , and 
the selection mode are used as inputs in a neural-network 
analysis trained to discriminate signal from background. 
The final selection requirement is a neural network out- 
put (NN) consistent with signal events. 

The signal B° — > t + t~ selection criteria for E res , NN 
and B pur are chosen to minimize the expected upper limit 
on B(B° — > t + t~). That optimization also rejects the 
signal selection modes t + t~ — > Ivv j pv and t + t~ — » 
-Kvjpv. After the full signal B° —> t + t~ selection, the 
combinatorial companion- B background is estimated and 



subtracted using ARGUS and Crystal Ball fits to the 
m ES distributions in simulation samples and data (Fig.[5J). 
From these fits we determine the signal efficiency (e S ig), 
the expected number of background events (N expecte d) , 
and the number of observed data events (N i, a ). Includ- 
ing systematic uncertainties described below, we obtain 
e sig = 0.043 ± 0.009, and N expected = 281 ± 48. We ex- 
tract from the fit 2V {, S = 263 ± 19 events in the data after 
the full selection. The central value of the B° — » t + t~ 
branching fraction is (—1.5 ± 4.4) x 10~ 3 . We find no 
evidence for signal events. Table ITU shows e S ig, N expecte d 
and N bs obtained from individual fits to specific signal 
selection modes. 

Systematic uncertainties on N expecte d and e S i g arise 
from several sources. The simulation statistical uncer- 
tainty for N expecte d (esig) is 10 events (11%). The system- 
atic uncertainties are estimated for cluster corrections to 
be 8 (3%), for particle identification corrections 10 (10%), 
and for tracking corrections 7 (3%). The m ES background 
subtraction fits after the full selection add a further sys- 
tematic uncertainty of 4 (2%). We estimate the system- 
atic uncertainty on N expecte d due to B decay modeling 
in EvtGen to be 10%. We estimate the systematic un- 
certainty due to our model of r decay by inserting dis- 
tributions obtained from the specialized r Monte Carlo 
TAUOLA to decay two r produced with the same 
helicity and the requisite momentum for a B° — > t + t~ 
decay. For each simulated event, the decay mode of each 
t is identified and the |p+|, |p_| and cos 9 values are re- 
placed with values sampled from distributions generated 
by TAUOLA for that mode. The relative e S i g variation 
between EvtGen and TAUOLA simulation is 2%. 

A final systematic uncertainty for both signal and 
background is assigned to the modeling of E res . The sim- 
ulation of background hits and hadronic interactions in 
the EMC does not perfectly model the data, and the dis- 
crepancy manifests itself in the E res distribution (Fig.|3Jl. 
This uncertainty is estimated from the difference between 
data and the simulation for a control process. The con- 
trol sample selection is identical to the B° — > t + t~ selec- 
tion except that events with an additional reconstructed 
Kg are selected and the K s daughters are removed from 
the event. For correct Kg reconstructions, this control 
sample models the K\ background while for Kg recon- 
structions from random combinations of tracks it models 
the backgrounds in which two oppositely charged parti- 
cles are lost due to the limited detector acceptance in 
the direction of the higher energy beam. The compo- 
sition of the background in the simulated control sam- 
ple agrees well with that of the simulated signal sam- 
ple. The control sample yields are 135 ± 14 events (data) 
and 125 ± 7 (simulation), for a relative discrepancy of 
(8 ± 13)%, consistent with zero. The systematic uncer- 
tainty due to modeling the residual energy in the EMC 
is taken to be the uncertainties in data and simulation 
yields added in quadrature, namely 13%. 
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FIG. 3: The E res distribution in the nominal sample (left) 
and the control sample (right) for data (dots), simulated back- 
ground (solid histogram), and simulated signal (dashed his- 
togram). The simulated signal distribution normalization is 
arbitrary. All requirements except those for E rea and NN 
are imposed. The events to the left of the vertical line are 
selected. 



Systematic uncertainties on the companion-S yield, 
expected background, and e S i g are folded into the upper 
limit calculation using the technique described in Ref. 
E3, giving 



B{B° 



T+T- 



) < 4.1 x 10" 



(4) 



at the 90% confidence level. The result constrains lep- 
toquark couplings as described in Ref. 2]. For example, 
the scalar SU(2) doublet leptoquark S1/2 can mediate 
B° — > t + t~. If no other leptoquark mediates the decay, 
the product of its coupling Af 3 (coupling righthanded b 
with r) with Ajj (coupling righthanded d with r) is 



A^A^ 3 < 1.4 x 10~ 2 



ms 1 



/2 



TOOGeVJ ' (5) 
at the 90% confidence level, where ms 1/2 is the S1/2 mass. 
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